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NONLINEAR EQUATIONS FOR FRACTIONAL 
LAPLACIANS II: EXISTENCE, UNIQUENESS, AND 
QUALITATIVE PROPERTIES OF SOLUTIONS 

Q . XAVIER CABRE AND YANNICK SIRE 

> 

q , Abstract. This paper, which is the follow-up to part I, concerns 

the equation (-A) s v + G'{v) = in 1", with s G (0,1), where 
(— A) s stands for the fractional Laplacian — the infinitesimal gen- 
erator of a Levy process. 

When n = 1, we prove that there exists a layer solution of the 
' equation (i.e., an increasing solution with limits ±1 at ±oo) if and 

, only if the potential G has only two absolute minima in [—1,1], 

located at ±1 and satisfying G'(— 1) = G'(l) = 0. Under the addi- 
tional hypothesis G"(— 1) > and G"(l) > 0, we also establish its 
uniqueness and asymptotic behavior at infinity. Furthermore, we 
provide with a concrete, almost explicit, example of layer solution. 

For n > 1, we prove some results related to the one-dimensional 
symmetry of certain solutions — in the spirit of a well-known con- 
1 jecture of De Giorgi for the standard Laplacian. 

ON 

t> 
q . 

1. Introduction 

This paper, which is a follow-up to our work [8], is devoted to study 
the nonlinear problem 

(-A) s v = f(v) inR", (1.1) 

where s G (0, 1) and 

(-A)M*) = C n , s P.V. / (1.2) 

is the fractional Laplacian. In the previous integral, P.V. stands for 
the Cauchy principal value and C n ^ s is a normalizing constant to guar- 
antee that the symbol of the resulting operator is |£| 2s ; see [S] for more 
details. As explained in section 3 of [S], this problem is equivalent to 
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the nonlinear boundary value problem 

div (y a Viz) = in ™ n+1 

: i + a )irT = /( M ) on arc 



du if \ ( 1 - 3 ) 

- — = flu) on dWl \ 
du a K ' + 



where n > 1, R+ +1 = {(a;,?/) G R n x R : y > 0} is a halfspace, 
<9R™ +1 = {y = 0}, u = u(x,y) is real valued, and 

<9« 1- on 

— hm ?/ a y it 



<9z/ a y^o ' 

is the generalized exterior normal derivative of it. Points in R™ are 
denoted by ). The parameter a belongs to (—1,1) and 

is related to the power of the fractional Laplacian (— A) s by 

a = 1 - 2s. 

Indeed, Caffarelli and Silvestre (see pUllE]) proved the following formula 
relating the fractional Laplacian (— A) s to the Dirichlet-to-Neumann 
operator: 

(711 

(-A) s {u(; 0)} = d s — in R n = dR n +\ (1.4) 

where d s is a positive constant depending only on s. 

The aim of the present paper is to study some special bounded so- 
lutions of (II. ip . The solutions we consider are the so-called layer solu- 
tions, i.e., those solutions which are monotone increasing, connecting 
— 1 to 1 at =Foo, in one of the x-variables. We focus on their existence, 
uniqueness, symmetry and variational properties, as well as their as- 
ymptotic behavior. 

In our previous paper [8], we proved a Modica-type estimate which 
allowed to derive a necessary condition on the nonlinearity / for the ex- 
istence of a layer solution in R. More precisely, we proved the following 
result. 

Theorem 1.1 (|8J). Let a G (—1,1) and f any C 1,7 (R) function, for 
some 7 > max(0, a). Let n = 1 and u be a layer solution of fll .3j) . that 
is, a bounded solution of (11.31) with n = 1 such that u x (-,0) > in R 
and u(x, 0) has limits ±1 as x — > ±oo. 

Then, for every x G R we have f Q °° t a \ Vu(x, t)\ 2 dt < oo and the 
Hamiltonian equality 

(1 + a) / -t a {u 2 x {x,t) -u 2 y {x,t)}dt = G{u{x,0)) -G(l). 
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Furthermore, for all y > and i6Kwe have 
ry i 

(l + o) / -t a {u 2 x {x,t)-u 2 y {x,t)}dt<G{u{x,0))-G(l). 
Jo ^ 

In the previous theorem, the last estimate is uniform as s tends 
to 1, i.e., as 1 + a tends to 0. This led in [8] to the convergence of 
layers, as s 1, to a layer of — v" = f(v) in R. In addition, using the 
Hamiltonian estimates of Theorem 11.11 we established the following 
necessary conditions for the existence of a layer in R. 

Theorem 1.2 (|8J). lei s G (0,1) and / any C 1 ' 7 (R) function, for 
some 7 > max(0, 1 — 2s). Assume that there exists a layer solution v 
of 

{-d xx ) s v = f(v) znR, (1.5) 
that is, v is a solution of (ll.5p satisfying 

v' > m R and lim v(x) = ±1. 
Then, we have 

G'(-l) = G'(l) = (1.6) 

and 

G > = G(l) in (-1,1). (1.7) 

In the present paper, we prove that the two necessary conditions 
in Theorem 11.21 are actually sufficient to ensure the existence of a 
layer solution in R. Under the additional hypothesis G"(—l) > and 
G"(l) > 0, we also prove the uniqueness (up to translations) of a layer 
solution in R and we establish its asymptotic behavior at infinity. 

To study the asymptotic behavior of the layer solution for a given 
nonlinearity /, it will be very useful to have the following almost ex- 
plicit example of layer solution for a particular nonlinearity. For every 
t > 0, a layer solution for some odd nonlinearity /* G C 1 ([— 1,1]) 
(see Theorem 13.11 below for more details) is provided by the following 
function: 

vl(x) = -l + 2 f Ps (t,x)dx=- f°° S ^le- tr2s dr, (1.8) 

J-oo 7T JO r 

where p s is the fundamental solution of the linear fractional heat equa- 
tion 

d t w + {-d xx ) s w = 0, t > 0, x G R. 
When s = 1/2, the particular layer solution above agrees with the 
explicit one used in [9], namely 

v i/2( x ) = ~ arctan — , with fy 2 (v) = —sm(irv). 



4 



XAVIER CABRE AND YANNICK SIRE 



In [9], J- Sola- Morales and one of the authors studied layer solutions of 



which corresponds to the case a = (that is s = 1/2) in ( 11. 3D . The 
goal of our paper is to generalize this study to any fractional power of 
the Laplacian between and 1. We will make a great use of the tools 
developed in [9]. 

The study of elliptic equations involving fractional powers of the 
Laplacian appears to be important in many physical situations in which 
one has to consider long-range or anomalous diffusions. From a proba- 
bilistic point of view, the fractional Laplacian appears as the infinites- 
imal generator of a Levy process (see the book of Bertoin [5]). In our 
case, as in [H] , we will concentrate on the problem (jl.3p and we will not 
consider probabilistic aspects. 

Problem (II .3p is clearly a degenerate elliptic problem concerning the 
weight y a . However, since a 6 (—1,1), the weight y a belongs to the 
Muckenhoupt class of Ai functions, i.e., it satisfies 



where w(x,y) = \y\ a and B denotes any ball in M n+1 . This fact allows 
to develop a regularity theory for weak solutions of (jl.3P ; see [8]. 

Another important property of the weight y a is that it just depends 
on the extension variable y and not on the tangential variable x. The 
equation is therefore invariant under translations in x, which allows the 
use of the sliding method to get uniqueness of layer solution in IR, as 
well as monotonicity of solutions with limits ±1 at ±oo. 

Remark 1.3. Another interesting problem is to consider the existence 
of monotone solutions of equation (II. ip connecting v(x%, x n ) at — oo 
to v(x 2 , x n ) at +oo where both v and v are solutions of (— A) s w = 
f(w) in M n_1 . We will not address this problem here, but we believe 
that the methods developed in the present paper (and in [9J [8j) allow 
to deal with this type of problem. 



Throughout the paper we will assume that the nonlinearity / is of 
class C 1,7 (1R) for some 7 > max(0, 1 — 2s). We will denote by G the 
associated potential, i.e., 





2. Results 
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The potential G is uniquely denned up to an additive constant. 

Let P s = P s (x,y) be the Poisson kernel associated to the operator 
L a = div (y a V), with a = 1 — 2s. We then have (see section 3 of [8]): 
for v a bounded Cf oc (R n ) function, v is a solution of (11. ip if and only if 

u(-,y) = P s (-,y) * v, 

a function having v as trace on <9M+ +1 , is a solution of (11.31) with / 
replaced by (1 + a)dj 1 f = 2(1 — s)dj 1 f. Recall that d s is the constant 
from (II. 4p . It turns out that 2(1 — s)d~ 1 has a positive limit as s | 
1. This is the reason why we wrote problem (ll.3p in [3] with the 
multiplicative constant 1 + a = 2(1 — s) in it; we wanted uniform 
estimates as s t 1. 

Let us recall some regularity results from [8]. The first one is Lemma 
4.4 of 0. 

Lemma 2.1 Q8J). Let f be a C 1,7 (1R) function with'j > max(0, 1 — 2s). 
Then, any bounded solution of 

(-A) s v = f{v) mR n 

is C 2,l3 (M. n ) for some < (5 < 1 depending only on s and 7. 

Furthermore, given sq > 1/2 there exists < < 1 depending only 
on n, s , and 7 — and hence independent of s — such that for every 
s > s , 

lkl|c7 2 ^(R») < C 

for some constant C depending only on n, sq, WfWc 1 ^, and ||u||i,<»(R n ) 
— and hence independent of s G (so, 1). 

In addition, the function defined by u(-,y) = P s (-,y) * v (where P s 
is the Poisson kernel associated to the operator L a ) satisfies for every 
s > s , 

IMIc^i^+T) + IIVx^H^^n+T) + \\D 2 X U\\ CI3 ^ TT) < C 

for some constant C independent of s G (so, 1), indeed depending only 
on the same quantities as the previous one. 

Following [9J, we introduce 

B+ = {(x,y) ER n+1 : y > 0,\(x,y)\ < R}, 

T° R = {(x,0) G dWl +1 : \x\ < R}, 

r+ = {(x,y)eR n+1 :y>0,\(x,y)\ = R}. 
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We consider the problem in a half-ball 

div (y a Vm) = in 5+ 

In the sequel we will denote by 

L a = div(y a V) 

the differential operator in ( |2.ip . Obviously, there is a natural notion 
of weak solution of (j2.1jl ; see [8]. 

We have the following regularity result (Lemma 4.5 of [8]). 

Lemma 2.2 (|8j). Let a G (-1, 1) and R > 0. Let ip e C a (T° R ) for 
some o G (0, 1) and w G L°°(B^ R ) R i/ 1 (i?^ R , y a ) fre a u>ea£; solution of 

L a u = m C 

— — = o? on 1 o R . 

Then, there exists (3 G (0, 1) depending only on n, a, and a, such 
that u G CP(B+) and y a u y G C P {B%). 

Furthermore, there exist constants C R and C R depending only on n, 
a, R, \\u\\ Loo{b + r) and also on |MU°°(r° K ) (for Or) and |Mlc-(r° fl ) (for 
C R ), such that 

~r+\ < G\ (2.2) 



it 



and 

IbXIUij) < c r- ( 2 - 3 ) 

Problem (12.11) has variational structure, with corresponding energy 
functional 

E b+W = I , -VW* + [ j^Giw), (2.4) 



r+ 2 / r o 1 + a 

ft 1 ° 



where G' = —f. This allows us to introduce some of the following 
notions. 

Definition 2.3. a) We say that u is a layer solution of (II. 3p if it is a 
bounded solution of (II. 3p . 

u xi > on dRl +1 , and (2.5) 
lim u(x, 0) = ±1 for every (x 2 , . . . , x n ) G W 1 ' 1 . (2.6) 

X\— >±oo 

Note that we will indifferently call layer solution a solution as above 
for problem (II. 3p or a solution v of equation (II. ip satisfying the same 
properties. 
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b) Assume that u is a C 13 function in M™ +1 for some (3 G (0,1), 
satisfying — 1 < u < 1 in IR" +1 and such that for all R > 0, 

! / B |Vu| 2 eL 1 (S+). 
We say that u is a local minimizer of problem (11. 3ft if 

E B +(u)<E B +(u + <tl>) 

for every R > and every C 1 function ^ in IR™ +1 with compact support 
in Bft U T° R and such that — l<u + tp<lin B^. To emphasize this 
last condition, in some occasions we will say that u is a local minimizer 
relative to perturbations in [—1,1]. 

c) We say that u is a stable solution of (11.31) if u is a bounded solution 
of (USD and if 

/ y a \ vei 2 -/ t^/'H^o (2-7) 

for every function £ G C 1 (1R™ +1 ) with compact support in IR™ +1 . 

It is clear that every local minimizer is a stable solution. At the 
same time, it is not difficult to prove that every layer solution u is also 
a stable solution — for this, one uses Lemma 16.11 below and the fact 
that positive solution of the linearized problem to (II. 3p . 

2.1. Layer solutions in R. The following result characterizes the non- 
linearities / for which problem (11. II) admits a layer solution in K. In 
addition, it contains a result on uniqueness of layer solutions. 

Theorem 2.4. Let f be any C 1,7 (M) function with 7 > max(0, 1 — 2s), 
where s G (0, 1). Let G' = —f. Then, there exists a solution v of 

(-d xx ) s v = f(v) inf 

such that v' > in M. and lim^-too v(x) = ±1 if and only if 

G'(-l) = G"(l) = and G> G{-\) = G{\) in (-1,1). (2.8) 

// m addition /'(— 1) < anc? /'(I) < 0, t/ien t/iis solution is unique 
up to translations. 

As a consequence, if f is odd and /'(±1) < 0, t/ien the solution is 
odd with respect to some point. That is, v{x + b) = —v(—x + b) for 
some b G R. 

Remark 2.5. The statement on uniqueness of layer solution also holds 
for any nonlinearity / of class C l {\— 1,1]) satisfying /'(— 1) < and 
f'iX) < 0. There is no need for /' to be C 7 ([— 1, 1]). Indeed, we will 
see that the proof follows that of [1] and thus only requires / to be 
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Lipschitz in [—1, 1] and nonincreasing in a neighborhood of —1 and of 
1. See also Lemma 5.2 of [9] where this more general assumption is 
presented. 

Note that a layer solution v = v(x), x G R, as in Theorem 12.41 
provides with a family of layer solutions of the same equation in R n . 
More precisely, for each direction e G R™, with |e| = 1 and e\ > 0, let 

v e (x 1 , . . . , x n ) := v({e, (x 1 , . . . , x n ))). 

Then, v e is a layer solution of 

(-A)V = f{v e ) inR n . (2.9) 

This fact is not immediate from the definition of the fractional Lapla- 
cian (11.21) through principal values in R and in R n — indeed, the inte- 
grals in R and in R n differ, but the normalizing constants C n>s in front 
make them agree. This fact — that v e solves (12. 9p — follows directly 
from the equivalence of problem (II. ip with the extension problem (11.31) 
and the fact that the constant d s in (11.41) is independent of the dimen- 
sion n. 

The equality G(— 1) = G(l) is equivalent to 

J f(s)ds = 0. 

Remark 2.6. Note that G may have one or several local minima in 
(—1,1) with higher energy than —1 and 1, and still satisfy condition 
(12. 8p . Such G will therefore admit a layer solution, hence a solution 
with limits —1 and 1 at infinity. Instead, such layer solution will not 
exist if G has a minimum at some point in (—1,1) with same height 
as —1 and 1. In particular, when G is periodic (as in the Peierls- 
Nabarro problem f(u) = sin(7ra), see [21]), the previous theorem proves 
that there exists no increasing solution connecting two non-consecutive 
absolute minima of G. 

In [12], with different techniques than ours it is proved that for po- 
tentials G with G'(-l) = G'(l) = and G > G(-l) = G(l) in (-1, 1), 
there exists a layer solution to equation (II. ip . We also refer to the 
interesting paper [12] where properties of ground state solutions are 
investigated. 

Our next result gives the asymptotic behavior of layer solutions. 

Theorem 2.7. Let f be any C 1,1 (M.) function with 7 > max(0, 1 — 2s), 
where s G (0, 1). Assume that /'(— 1) < ; f'{l) < ; and that v is a 
layer solution of 

{-d xx ) s v = f(v) m R. 
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Then, there exist constants < c < C such that 

c\x\- l ~ 2s < v'(x) < C\x\- X ' 2s for \x\ > 1. (2.10) 
As a consequence, for other constants < c < C , 

cx~ 2s < 1 - v(x) < Cx~ 2s for x>\ (2.11) 

and 

c\x\~ 2s < 1 + v(x) < C\x\- 2s for x < -1. (2.12) 

To prove the above theorem for a given nonlinearity /, the following 
almost explicit layer solution (we emphasize that it is a layer solu- 
tion for another nonlinearity) will be very useful. More properties and 
remarks on this concrete layers will be given in section 3. 

Theorem 2.8. Let s G (0, 1). For every t > 0, the C°°(R) function 
vl(x) := -1 + 2 f p t (t,x) dx=- [°° !^M e -^ dr 

J-oo * J0 r 

is the layer solution in K of (11. ip for a nonlinearity /* G C 1 ([— 1, 1]) 
which is odd and satisfies fl(0) = = 0, /* > in (0,1), and 

(/i)'(±i) = 1 1- 

In the theorem, since /* G ^([-l, 1]) and (/ s *)'(±l) < 0, Theo- 
rem 12.41 and Remark 12.51 guarantee that its corresponding layer v * is 
unique up to translations. 

As we will see in Theorem 12.111 below, every layer solution is a lo- 
cal minimizer and, in particular, a stable solution. This holds in any 
dimension and for any nonlinearity. Our next result states that the 
converse is also true in dimension one and under certain hypothesis on 
the nonlinearity. That is, under various assumptions on G, we prove 
that, for n = 1, local minimizers, solutions with limits (not monotone 
a priori), or stable solutions are indeed layer solutions. 

Theorem 2.9. Let f be any C ,1,7 (R) function, with 7 > max(0, 1 — 2s). 
Let n = 1 and u be a function such that 

\u\ < 1 in R\. 

a) Assume that G > G(—l) = G(l) in (—1, 1), and that u is a local 
minimizer of problem (11.31) relative to perturbations in [—1,1]. Then, 
either u = u(x,y) or u* = u*(x,y) := u(—x,y) is a layer solution 

of (US). 

b) Assume G"(—l) > 0, G"(l) > 0, and that u is a solution of (11.31) 
with 

lim u(x, 0) = ±1. 

x—>±oo 
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Then, u is a layer solution of ( II. 3p . 
c) Assume that G satisfies: 

if - 1 < L~ < L + < 1, G'(L-) = G'(L + ) = 0, (2.13) 
and G > G(L-) = G(L + ) in (L",L + ), (2.14) 
then IT = -1 and L+ = 1. (2.15) 

Let u be a nonconstant stable solution of ( 11.3j) . Then, either u = 
u(x,y) or u* —u*(x,y) :— u(—x,y) is a layer solution of (11. 3p . 

Remark 2.10. Notice that the hypothesis ( I2.13p - ()2.15p on G in part c) 
of the theorem is necessary to guarantee that u connects ±1. Indeed, 
assume that — 1 < L~ < L + < 1 are four critical points of G with 
G > G(-l) = G(l) in (-1,1) and with G > G{L~) = G{L + ) in 
(L~,L + ). Assume also that 

G{±1) < G(L ± ). 

Then, by our existence result (Theorem 12. 4p applied twice — in (—1, 1) 
and also in (L~, L + ) after rescaling it — , we have that (— d xx ) s v = f(v) 
in M admits two different increasing solutions: one connecting L ± at 
infinity, and another connecting ±1. 

Instead, as pointed out in Remark 12.61 if G > G(±l) = G^L^ 
in (—1,1), then there is no increasing solution connecting ±1, as a 
consequence of our Modica estimate, which gives (II. 7p . 

Note that an identically constant function u = s is a stable solution 
of dUD if and only if G'(s) = and G"(s) > 0. This follows easily 
from the definition (12.71) of stability. Therefore, regarding part c) of 
the previous theorem, a way to guarantee that a stable solution u is 
nonconstant is that u = s G (—1, 1) at some point and that either 
G'(s) ^ or G"(s) < 0. 

2.2. Stability, local minimality, and symmetry of solutions. 

The following result states that every layer solution in R^ +1 is a lo- 
cal minimizer. This result is true in every dimension n. 

Theorem 2.11. Let f be any C 1 ' 7 (R) function and'y > max(0,l — 2s), 
where s G (0, 1). Assume that problem (II. 3p admits a layer solution u 
in M™ +1 with n > 1 . Then : 

a) u is a local minimizer of problem (11.31) . 

b) The potential G satisfies 

G'(-l) = G"(l) = and G > G(-l) = G(l) in (-1, 1). (2.16) 



The strict inequality G > G{— 1) = G{1) in (12.161) is known to hold 
when n = 1 or, as a consequence, when w(-,0) is a one- dimensional 
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solution in R n . We established this in [8] (it is one of the implications 
in Theorem 12.41 above). The strict inequality G > G(±l) also holds 
when n = 2 (as a consequence of Theorem 12. 121 below) and when n = 3 
and s > 1/2 (as a consequence of a result from [7]). It remains an open 
question in the rest of cases. 

For n = 2, we prove that bounded stable solutions u (and hence also 
local minimizers and layer solutions) are functions of only two vari- 
ables: y and a linear combination of X\ and x 2 - This statement on the 
ID symmetry of u(-, 0) is closely related to a conjecture of De Giorgi on 
ID symmetry for interior reactions, proved in |13[ El 12] in low dimen- 
sions and partially settled by Savin [T8] up to dimension 8. We also 
refer the reader to |20] where some rigidity properties of boundary 
reactions have been established through a more geometric approach. 
Particularly, in [19], the following symmetry result in dimension n = 2 
is proved by using a completely different approach than the one used in 
the present paper, relying on a weighted Poincare inequality (see also 

EH)- 

Theorem 2.12. Let f be anyC ln (R) function and'-f > max(0, 1 — 2s), 
where s G (0, 1). Let v be a bounded solution of 

(-A) s v = f(v) in R 2 . 

Assume furthermore that its extension u is stable. 

Then, v is a function of one variable. More precisely, 

v(xi,x 2 ) = v (cos(6 l )xi + sin(6 , )x 2 ) in R 2 

for some angle 9 and some solution Vq of the one- dimensional problem 
with same nonlinearity f , and with either v' > everywhere or vq 
identically constant. 

For n = 3 and s > 1/2, this ID symmetry result has been proved 
by E. Cinti and one of the authors in [61 Ej. It remains open for n = 3 
and s < 1/2, and also for n > 4. 

A simpler task than the study of all stable solutions consists of study- 
ing solutions u of f ll.3p with \u\ < 1 and satisfying the limits f !2.6p 
uniformly in (x2, ■ ■ ■ ,x n ) G M n_1 . Under hypothesis /'(— 1) < and 
/'(— 1) < 0, it is possible to establish in every dimension n that these 
solutions depend only on the y and x\ variables, and are monotone 
in x\. Here, by the uniform limits hypothesis, the x- variable in which 
the solution finally depends on is known a priori — in contrast with the 
variable of dependence in Theorem 12.121 For the standard Laplacian 
this result was first established in [3] using the sliding method. We will 
not provide the proof of the result because it is completely analogue 
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to the one in [9]. Since our operator L a is invariant under translations 
in x, one can perform the sliding method together with the maximum 
principles proved in [8]. 

Theorem 12. 12l is a partial converse in dimension two of Theorem 12. Ill 
a), in the sense that it establishes the monotonicity of stable solutions 
and in particular, of local minimizers. The remaining property for being 
a layer solution (i.e., having limits ±1 at infinity) requires additional 
hypotheses on G, as in Theorem 12.91 

2.3. Outline of the paper. In section 3 we construct an almost ex- 
plicit layer solution (Theorem 12. 8p and we use it to establish the as- 
ymptotic behavior of any layer solution in R as stated in Theorem 12.71 
In section 4 we prove the existence of minimizers to mixed Dirichlet- 
Neumann problems in bounded domains of — a result needed in 
subsequent sections. In section 5 we prove the local minimality of layer 
solutions in any dimension and the necessary conditions on G for such a 
layer in R n to exist, Theorem 12. Ill The ID symmetry result for stable 
solutions in M 2 , Theorem 12. 121 is established in section 6. Finally, sec- 
tion 7 concerns layers in R and establishes the existence Theorem 12.41 
and the classification result Theorem 12.91 



3. An example of layer solution. Asymptotic properties 

of layer solutions 

In this section we provide with an example of layer solution based on 
the fractional heat equation. From it, we get the asymptotic behavior of 
layers for all other nonlinearities. Let us first explain how the concrete 
layer is found. 

The starting point is the fractional heat equation, 

d t w + {-d xx ) s w = 0, t > 0, x e R, (3.1) 

which is known to have a fundamental solution of the form 

Ps {t,x) =r&q t (t-&x) > (3.2) 

for a;6l,t>0. Being the fundamental solution, p s has total integral 
in x equal to 1, i.e., 

p s (t, x) dx = 1 for all t > 0. (3.3) 

To compute p s , one takes the Fourier transform of (13.11) to obtain 

dtPs + \Z\ 2s Ps = 0, 
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where p s = p s (t, £) is the Fourier transform in x of p s (t, x). Thus, since 
p s (0, ■) is the Delta at zero and hence p s (0, •) = 1, we deduce 

p s (t,0=exp{-m 2s }. 

From this, by the inversion formula for the Fourier transform, we find 

1 f°° 2 

p s (t,x) = — cos(xr)e~ tr s dr. (3.4) 
ir Jo 

It follows that the function 

/X r-X 
p s (t,x)dx = 2 p s (t,x)dx (3.5) 
■oo J0 

is increasing and has limits ±1 at ±oo. The concrete expression (13. 6 p 
below for u* is obtained by interchanging the order of the two integrals 
when using (13 .4p to compute the primitive of p s . 

That v\ is a layer solution is stated in the next theorem, which con- 
tains all statements in Theorem 12.81 and also the asymptotic behavior 
of vl, among other facts. The proof of the theorem is given at the end 
of this section. 

Theorem 3.1. Let s G (0, 1). For every t > 0, the C°°(M) function 

PX PX 

vl(x) := — 1 + 2 / p s (t,x)dx = 2 j p s (t,x)dx 

-e~ tr ^ dr 



2 f°° sin(xr) _.-zs 



7T In r 



szgn(x) - [°° ^i e -*(*/M) 2s dz (3.6) 



7T 







is the layer solution in K of (11. ip for a nonlinearity /* G C 1 ([— 1, 1]) 
which is odd and twice differentiable in [—1, 1] and which satisfies 



fM = f s (l) = 0, /*>0ro(0,l), {fl)\±l) = -- t 



and 



. < z/0<s<l/2, 
4 7T cos(tts) r(4s) I 

(/J( 1 ) = -T— 7 — r / r / u 2 S = if s = 1/2, (3.7) 
t sm(7rs) (r(2s))^ ' 



>0 i/l/2<s<l. 

In addition, the following limits exist: 

4s 

lim |x| 1+2s (<9 x ^)(x) = t— sin(7rs)r(2s) > (3.8) 

|x|->DO 7T 
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and, as a consequence, 

2 

lim |x| 2 >*(x) =F 1| = t-sin(7rs)r(2s) > 0. 

x— »±oo 7T 

Remark 3.2. As stated in the theorem, we have /* G C 1 ([— 1, 1]) and 
(/*)'(±1) < for every s G (0,1). In particular, by Theorem 12.41 and 
Remark I2.5[ its corresponding layer v\ is unique up to translations. 

When s = 1/2, the particular layer above agrees with the explicit 
one used in [9], namely 

v i/2( x ) = ~ arctan — , with fi/ 2 ( v ) = — sin(7rt>). 

This can be easily seen computing fl3.4[) explicitly when s = 1/2, using 
integration by parts, to obtain 

We may try to see which function we obtain in the above formulas 
setting s = 1. In this case, (\3A\i can be checked to be equal to a 
Gaussian and thus v\, two times its primitive, is the error function 
erf(x) — up to a scaling constant. Its derivative is therefore e~ cx , 
which does not have the correct decay e~ cx at +oo for the derivative 
v' of a layer solution to — v" = f{v). This is due to the fact that the 
limit as s — >• 1 of /* will not be a C 1 ([— 1, 1]) nonlinearity at the value 
1 — even if they all satisfy (/*)'(1) = —1/t. The reason is that their 
second derivatives at 1, (/*)"(1), blow-up as s —> 1 as shown by (13. 7p . 

Note also that (13.71) shows that, when 1/2 < s < 1, the nonlinearity 
/* is positive but not concave in (0, 1). 

The following immediate consequence of Theorem 13.11 will give the 
asymptotic behavior of layer solutions for any nonlinearity /. 

Corollary 3.3. Let s G (0, 1) and t > be a constant. Then, the 
function 

(p f = d x v\ > 0, 
where v\ is the explicit layer of Theorem \3.1\ satisfies 

(-d xx ) s ip\x) + 2t~ l ip\x) > for \x\ large enough, (3.9) 

(—d xx ) s <p t (x) + 2~ 1 r l ip\x) < for \x\ large enough, (3.10) 
and that the following limit exists and is positive: 

lim |x| 1+2 V(x) G (0,+oo). (3.11) 

|x|— >oo 
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Proof. Clearly = d x v\ > satisfies the linearized equation 

(-^)V-(/ s t )K(x))^ = inR. 

Using that y?' > 0, v l s has limits ±1 at ±oo, /* is C 1 ([— 1, 1]), and that 
(/*)'(±1) = -1/t, both §3M) and fkTTOj) follow. The statement (SHI) 
follows from (EES)) . □ 

With this corollary at hand, we can now prove the asymptotics of 
any layer. 

Proof of Theorem 12.71 The proof uses Corollary 13.31 above and a very 
easy maximum principle, Lemma 4.13 and Remark 4.14 of Its 
statement in dimension one is the following. 

Let w e C^ C (R) be a continuous function in R such that w(x) — » 
as \x\ —> oo and 

{-d xx ) s w + d(x)w > inR (3.12) 

for some bounded function d. Assume also that, for some nonempty 
closed set H c R, one has w > in H and that d is continuous and 
nonnegative in R \ H. Then, w > in R. 

Let now / and v be a nonlinearity and a layer as in Theorem 12.71 
We then have 

(-d xx ) s v'-f'(v)v' = Q inR. (3.13) 
To prove the upper bound for v' in (|2.10|) . we take t large enough such 
that 2t _1 < min{— /'( — 1), — /'(l)}. Then, for any positive constant 
C > 0, 

w := CV* - v' 

satisfies, by (I3.9P and f)3.13p . (— d xx ) s w+2t~ 1 w > for \x\ large enough, 
say for x in the complement of a compact interval H. Next, take the 
constant C > so that w > 1 in the compact set H, and define now 
d in if so that (— d xx ) s w + cfoy = in H — recall that w > 1 in H and 
hence d is well defined and bounded in H . We take d = 2t~ l in R \ H. 
Thus, (I3.12p is satisfied and, since w — > at infinity, the maximum 
principle above leads to w > in R. This is the desired upper bound 
for v' in (12.101) . since ip l satisfies (13.111) . 

To prove the lower bound for v' in (I2.1UI) . we proceed in the same 
way but replacing the roles of v' and </?*. For this, we now take t > 
small enough such that max{— /'(— 1), — /'(l)} < 2 _1 t _1 . Thus, w : = 
Cf' — (p* satisfies (— d xx ) s w + 2~ 1 t~ 1 w > for |x| large enough. One 
proceeds exactly as before to obtain w > in R for C large enough, 
which is the desired lower bound for v' in (12.101) . □ 

It remains to establish Theorem 13.11 For this, we use the following 
well-known technical lemma due to G. Polya [17J, 1923. We prove it 
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here for completeness; in fact, the proof as explained in [T7] only works 
for s < 1/2. For s > 1/2, we follow the proof given in 



Lemma 3.4. For n > and s G (0, 1), we have 

POD 

lim / sm(z)z ns ~ 1 e~ i ~ z/x)2a dz = siu.{ksti /2)T(ks). 

x^+od J Q 

Proof. For every x > 0, we have 



oo 



sm(z)z KS 1 e (z/:r) 8 dz = Im / e iz z KS 1 e lyZlx) 8 dz 
o Jo 



where 

Let us also denote 



poo 

= Im / h x (z) dz, 
Jo 

h x (z) := ^«-i e «-(*/*) 2 *. 



For < 9 < 7r/2, let 70 be the half- line from the origin making an 
angle 9 with the positive x-axis. We will next see that, for certain 
angles 9, Im f h x (z) dz are all equal and independent of those 9. For 
this, given two angles < Q\ < 9 2 < 7r/2 and R > 0, we integrate 
counter-clockwise on the contour given by the segments of length R 
starting from on 7^ and on 7# 2 , and by the arc e of radius R 
with center at the origin and joining the two end points of the previous 
segments. We also need to remove a neighborhood of zero and add a 
small arc with center at the origin connecting the two half-lines. The 
integrals of h x and of in this small arc will tend to zero as the radius 
tends to zero, since l^l- 2 )! + |^oo(^)| < C , |,z| K ' s ~ 1 near the origin. 

The key point is to make sure that the integral of h x , and later of 
/ioo, on the arc Tg 62 of radius R tends to zero as R — > 00 if we choose 
the angles < Q\ < 9 2 < 71/ 2 correctly. Note that if z G C belongs to 
such an arc, then z belongs to the sector 

S du e 2 := {z G C : d 1 < Arg(z) < 9 2 }. 

To guarantee the convergence to zero of the integral on the arc, note 
that 

= |z| KS - 1 exp{-Im(z) - x~ 2s Re(z 2s )} (3.14) 

and 

\hoo(z)\ = {zl™- 1 exp{-Im(z)} (3.15) 

for all z G C in the first quadrant. 
We need to distinguish two cases. 
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Case 1. Suppose that s < 1/2. In this case we take #i = and 
9 2 = 7r/2. Then, if z lies in the sector (the first quadrant), then 

z 2s is also in the first quadrant, since 2s < 1. Thus, the real and 
imaginary parts appearing in (I3.14p are both nonnegative, and at least 
one of them positive up to the boundary of the quadrant. Thus, by 
(I3.14p . \h x \ —> exponentially fast — as exp{— c(x)|2;| 2 ' s } — uniformly 
in all the quadrant. Hence, the integral on the arc ^Q n / 2 tends to zero 
as R — > oo. We deduce 

poo p r 

I sm{z)z KS - l e^ zlx)23 dz = lm / h x (z) dz = Im / h x (z) dz 

JO J 70 J It, /2 

= lmL iKSn/2 [ y KS - l e- y - l2s{y/x)2s dyX . 



Note that the function in the last integral is integrable since 

\ e -y-i 2s (y/x) 23 | _ | e -J/-(cos(s7r)+isin(s7r))(i//a;) 2s | __ e -y-cos{sn)(y / 'x) 2s < g 

due to s < 1/2. Thus, the limit as x — > +oo exists and is equal to 



/■OO f POO \ 

lim / sin(z)z KS - 1 e- (z/x)2s dz = lm\e iKSn/2 y KS - l e~ y dy\ 

x-++aoj { J Q J 

poo 

= sin(/ts7r/2) / y KS - x e- y dy 
Jo 

= sm(KSTt/2)T(Ks), (3.16) 

as claimed. 

Case 2. Suppose now that 1/2 < s < 1. In this case (I3.14p does not 
tend to zero at infinity in all the first quadrant, since 2s > 1 and thus 
Ke(z 2s ) becomes negative somewhere in the quadrant. Here, we need 
to take 

0i = and 6 2 = — . 

4s 

Now, in the sector So ;7r /(4 S ), the real and imaginary parts appearing in 
(13 . 141) are both nonnegative, and at least one of them positive up to 
the boundary of the sector. Thus, as before, we now deduce 

POO P P 

/ sin^^-V^ 28 dz = Im / h x {z) dz = \m I h x {z) dz. 

JO "'70 Jl7T/(4s) 

Note that in the last integral on j^/Us), we have 

1^(^)1 = |z| KS - 1 exp{-Im(z) - x- 2s Re(z 2s )} 
= |z| KS " 1 exp{-Im(z)} = |/ioo(^)| 
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for z G 77r/(4s)- Besides, by the last expression, is integrable on 
7tt/(4s)- Thus, by dominated convergence, we have 

/*oo r 

lim / sin(z)z KS - l e^ z/x)23 dz = Im / h^z) dz. (3.17) 

Finally, for this last integral we work on the sector SV/^s),^- By 
(I3.15p . h^z) tends to zero exponentially fast and uniformly as \z\ — > 00 
on the sector. Thus, 

Im / h co (z) dz = Im / h oa (z) dz 

Jin/Us) JItt/2 

= Im|e^ /2 ^ y KS ~ l e- y dy^. 

Recalling (13.171) . one concludes as in (I3.16p . □ 

Finally, we can prove our results on the explicit layer. 

Proof of Theorem \3.1\ Let v\ be defined by (13. 6p . It is clear that 

v t s (x)=v 1 s (t- l/i2s) x). (3.18) 

Hence, by the definition (11.21) of the fractional Laplacian, we have 
(— d xx ) s vl(x) = ^(—dxxYvKt' 1 ^ 23 ^). Thus, having proved all the 
statements for vl, they will also hold for v\ with nonlinearity f l s {v) = 

Hence, we may take t — 1. To simplify notation, we denote 



v 



and / := fl 



From v'(x) = 2p s (l, x) and expression (13. 4p . it is clear that v G C c 
By expression (I3.5p . we have v(— 00) = — 1. Since v' = 2q s = 2p s (l, ■) > 
0, v is increasing. 

The fact that v(+oo) = 1 is a consequence of (13.31) . J R p s (l, y) dy = 1. 
It also follows from expression (13. 6p and the well-known fact that, in 
principal value sense, J °° sin(z)z _1 dz = n/2. This can also be proved 
adding a factor z KS in the integral (13. 6p , using then Lemma 13.41 and 
noting that sin(/?S7r / '2)T (ks) = sin(ftS7r /2)(k,s)~ 1 T(hs + 1) — > n/2 as 

We now prove that there exists a function / such that 

{-d xx ) s v = f(v) inK. 

For this, we use the expression (13. 2p and that p s solves the fractional 
heat equation ( 13. ip . Because of the commutation of the derivative with 
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the fractional Laplacian, we deduce 

{(-cUVK(x) = (-cyV(rr) = 2{-d xx ) s q s {x) = -2d tPs {l,x) 

= -{q a (x) + xq' s (x)}. 
s 

Therefore, integrating by parts, 

1 f x , 1 1 

(-d xx ) s v (x) = - {q s (z) + zq' s (z)} dz = -xq s (x) = —xv'(x). 

S J-oo s 2s 

Since v' > 0, the C°° function v = v(x) is invertible on IR, with 
inverse x — x(v), a C°° function on the open interval (—1,1). We now 
set 

f(y) := ±x(v)v'(x(y)), (3.19) 

so that our semilinear fractional equation is satisfied. We know that 
/ G C°°{— 1, 1). Also, since v is an odd function, its inverse x is also 
odd and therefore / is odd, by H3. 191) . This expression also gives that 
/>0in(0,l). 

It remains to verify that / G C 1 ([— 1, 1]) once we set /(±1) = and 
/'(±1) = —1, and that / is twice differentiable in [—1, 1] and having 
values for /"(±1) given by f !3.Tj) with t — 1. It also remains to establish 
the asymptotic behavior of v' . 

For all this, using H3.4[) we compute 

{tt/2)v'(x) = Tiq s {x) = - / cos(^)e- ( ^ )2s dz 

X Jo 

1 f°° 

= - / {sm(z)}'e- izlx)2s dz 
x Jo 

POO 

= 2.SX- 1 - 23 / sin(z)z 2s - 1 e- (z/x)2fl dz, (3.20) 
Jo 

by integration by parts. Hence using Lemma [33] with k = 2, we deduce 

4s 

lim x 1+2s v'(x) = lim 2x 1+2s q s (x) = — sin(7rs)r(2s), (3.21) 

x— ¥+oo x— >+oo 7T 

as claimed in ( 13. 8ft — for other values of t, simply use ( 13. 18ft . In par- 
ticular, lim x _j. +00 xv '(x) = and thus, by ( 13.191) . / is continuous on 
[—1, 1] defining /(±1) = 0. In addition, we also deduce 

f°° 2 

1 - v (x) = / v'(y) dy = - sin(7rs)r(2s)x" 2s + o(x" 2s ) (3.22) 

Jx 7T 

as x — > +oo. 
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Next, we differentiate (I3.19p . that is, f{y(x)) = (2s)~ 1 xv'(x) 
(2s)~ 1 x2q s (x), to obtain 

f(v)v' = ^{v' + x(v)2q' s (x(v))} 



and hence 



/» = l{l + lW ^M)j. (3.23) 
2s [ q s {x{v)) ' 



Thus, using (13. 4 1) we compute 

/•oo 

7rxg^(x) = — / sin(xr)e _r s c/r 



/>oo 

= -aT 1 / zsin^e - ^^ cfe 

POO 

= -x 1 {sintz) - zcos(z)}'e~( z/x ^ 2s dz 
Jo 

POO 

= -2.SX- 1 - 28 {sin(z) - zcos{z)}z 2s - l e- {z/x) ' 2s dz. 
Jo 

We also compute vr{(l + 2s)q s + xq' s } by adding ( 13 . 20 j) (multiplied by 
1 + 2s) to the previous expression. Integrating by parts, and at the end 
invoking Lemma [3.41 with k = 4s, we obtain 



7i{{l + 2s)q s + xq' s } = 

POO 

= 2SX- 1 - 23 {2ssm(z) + zcos(z)}z 2s ~ 1 e^ z/x)2s dz 
Jo 

poo 

= 2sx^- 2s / {sin(^ 2s }V (2/x)2s dz 
Jo 

POO 

= (2s) 2 x- l - 4s / sin^z^-V^ 2 * dz 
Jo 

= x- 1 - 4s {4s 2 sin(2 7 rs)r(4s) + o(l)} 

= a;- 1 " 4s {8s 2 sin(7rs)cos(7rs)r(4s) +o(l)} (3.24) 
as x +oo. 

Therefore, from (TP3l . (l3~24l) . and <K21\\ . one has 

f ivix)) = _! + _L + + = -1 + 0(«-*). (3.25) 
2s q s 

Thus, setting /'(±1) = —1 and using that /' is even, we have that / 
is differentiable at ±1. 
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Finally, using (13T251) . ( I3~24l) . (13T22|) . and (13T2T1) . we have 

f (x) — 1 v(x) — 1 

1 (l + 2s)<? s + x^ cos(vrs)r(4s) 

— > —71- 



2s (v-l)q s sin(7rs)(r(2s)) 2 

as x — > +oo. This establishes that / G C x (\—1, 1]) and also that / is 
twice different iable in all of [—1, 1] with 

,///,-, x cos(7r£) T(4g) 

/ (±1) = ^7T— — - (3.26) 
sm(7rs) (r(2s)) 2 

The proof is now complete. □ 

4. MlNIMIZERS OF THE DlRICHLET-NEUMANN PROBLEM IN 

BOUNDED DOMAINS 

In this section, we concentrate on the existence of absolute mini- 
mizers of the functional E^{u) on bounded domains Q. This is an 
important step since, as in [9], the existence theory of layer solutions 
goes through a localization argument in half-balls of M" +1 . 

Let fl C be a bounded Lipschitz domain. We define the fol- 

lowing subsets of dQ: 

d°n = {(x, 0) G dWl +1 : B+(x, 0) C ft for some e > 0} and (4.1) 



d + n = dnnwf 1 . (4.2) 

Let H 1 (Q,y a ) denote the weighted Sobolev space 

H\Q,y a ) = {u-.Q^R : y a {u 2 + \Vu\ 2 ) G L\Q)} 

endowed with its usual norm. 

Let u G C"(p) R H 1 (Q,y a ) be a given function with |u| < 1, where 
/3 G (0, 1). We consider the energy functional 

En(v) = [ y -^\Vvf + [ j^-G(v) (4.3) 

in the class 

Cu,o(ft) = {f G H l (tt,y a ) :-1<v<1 a.e. in ft and u = u on <9 + ft}, 

which contains u and thus is nonempty. 

The set C u<a (VL) is a closed convex subset of the affine space 

H u , a (n) = {«£ H l {yi iy a ) : v = u on <9+ft}, (4.4) 

where the last condition should be understood as that v — u vanishes 
on <9 + ft in the weak sense. 
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Lemma 4.1. Let n > 1 and C 6e a bounded Lipschitz domain. 

Let u G (7^(0) Pi H 1 {yL 1 y a ) be a given function with \u\ < 1, where 
G (0, 1). Assume that 

f(l) < < /(-l). (4.5) 

T/ien, the functional Eci admits an absolute minimizer w inC u ^ a {Vt). 
In particular, w is a weak solution of 

L a w = in n 

(l + a)^ = /H ond°Q (4.6) 
w = u on d + £l. 

Moreover, w is a stable solution of (14. 6p . in the sense that 

y a m 2 - [ -i-/(«,)e 2 > o (4.7) 

/or every £ G H l (Q,y a ) such that £ = on d + Q in the weak sense. 

Hypothesis ( 14.51) states simply that —1 and 1 are a subsolution and 
a supersolution, respectively, of (14. 6ft . 

Proof of Lemma 14. 11 As in [9 J , it is useful to consider the following 
continuous extension f of f outside [—1, 1]: 



/(*) = 

Let 

G(s) = - I f, 
and consider the new functional 

Ea(v)= [ 7^1* + / T— G(G 




in the affine space H Ui0i (fl) defined by (I4.4p . 

Note that G = G in [—1, 1], up to an additive constant. Therefore, 
any minimizer w of Eq in H u ^ a (Q) such that —1 < w < 1 is also a 
minimizer of Eq in C u ^ a (Q). 

To show that Eq admits a minimizer in H u a (Q), we use a standard 
compactness argument. Indeed, let v G H Uta (Q). Since v — u = 
on d + Q, we can extend v — u to be identically in \ fi, and we 
have v — u G H 1 (Rl +1 ,y a ). By Nekvinda's result [15], the trace space 
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of H 1 (Rl +1 ,y a ) is the Gagliardo space W^' 2 {R n ) = H s (R n ). The 
Sobolev embedding (see [T]) 

H'(R n ) m> L^(M. n ) 

(or into any L P (M") if n = 1 < 2s) gives then the compactness of the 
inclusion 

H V)a {n) CC L 2 (d°Q). 

Now, since H u>a (Q) C L 2 (<9°fi) and G has linear growth at infinity, it 
follows that is well defined, bounded below, and coercive in H Uia (fl). 
Hence, using the compactness of the inclusion H U)a (Q) CC L 2 (d°Q), 
taking a minimizing sequence in H Uya {Vl) and a subsequence convergent 
in L 2 (d°Q), we conclude that Eq admits an absolute minimizer w in 

Since / is a continuous function, E is a C 1 functional in H Uta (Jl). 
Making first at second order variations of E at the minimum w, we 
obtain that w is a weak solution of (14.61) which satisfies (14.71) . with / 
and /' replaced by / and /', respectively, in both (14.61) and (14. 7p . 

Therefore, it only remains to show that the minimizer w satisfies 

— 1 < w < 1 a.e. in Q. 



We use that —1 and 1 are, respectively, a subsolution and a superso- 
lution of (I4.6p . due to hypothesis (14.51) . We proceed as follows. We 
use that the first variation of Eq at w in the direction (w — 1) + (the 
positive part of w — 1), is zero. Since \w\ = \u\ < 1 on d + Q and hence 
(w — 1) + vanishes on d + Q, we have that w + e(w — 1) + G H ua (Q) for 
every s. We deduce 

0= I y a VwV{w- l) + - [ f(w)(w-l) + 
Jq Jd a n 

y°\V(w-l) + \ 2 - [ /(1)(^-1) + 

Qn{w>l} Jd°Qn{w>l} 



> / y a \V(w-l) 



1 2 
I i 



n 



where we have used that f(s) = /(l) for s > 1, and that /(l) < 
by assumption. We conclude that (w — 1) + is constant, and hence 
identically zero. Therefore, w < 1 a.e. The inequality w > — 1 is 
proved in the same way, now using /(— 1) > 0. □ 
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5. Local minimality of layers and consequences. 
Proof of Theorem 12.111 

The fact that for reactions in the interior (that is, s — 1 in our 
equation), layer solutions in R™ are necessarily local minimizers was 
found by Alberti, Ambrosio, and one of the authors in [2]. For the 
fractional case, this is the statement in Theorem 12.111 a) above. The 
proof in [2] also works in the fractional case, working with the extension 
problem. It uses two ingredients: the existence result from the previous 
section (Lemma l4.ll) and the following uniqueness result in the presence 
of a layer. 

Lemma 5.1. Assume that problem (11. 3p admits a layer solution u. 
Then, for every R > 0, u is the unique weak solution of the problem 

' L a w = zn5+cM" +1 
-1<w<1 in 

dw_ , ^ 
dv a 



: i + «)7T37 = /H onT° R 



= u on F^. 

Proof. We refer the reader to the proof of Lemma 3.1 in [9 J since the 
proof is identical in our case. Indeed, since the operator L a is invariant 
under translations in x, this allows to use the sliding method as in 
Lemma 3.1 of [9] to get the uniqueness. The only other important 
ingredient in the proof is the Hopf boundary lemma; in our present 
context it can be found in Proposition 4.11 and Corollary 4.12 of [8]. □ 

Part b) of Theorem 12 . 1 1 1 will follow from the following proposition. It 
will be useful also in other future arguments. Notice that the result for 
n — 1 follows from our Modica estimate, Theorem 2.3 of [8] (rewritten 
in Theorem 11.11 of the present paper). Instead, the following proof 
also works in higher dimensions but only gives G > G(L~) = G(L + ) 
in [—1, 1] — in contrast with the strict inequality G > £?(— 1) = G(l) 
obtained in dimension one from the Modica estimate when = ±1. 

Proposition 5.2. Let u be a solution of (11. 3p such that \u\ < 1, and 



lim u(x, 0) = L for every (x 2 , • • • , %n) £ 

X\— >±oo 



on— 1 



for some constants L and L + {that could be equal). Assume that u is 
a local minimizer relative to perturbations in [—1, 1]. Then, 

G > G(L~) = G(L + ) in [-1,1]. 

Proof. It suffices to show that G > G(L~) and G > G(L + ) in [-1,1]. 
It then follows that G{L~) = G{L + ). By symmetry, it is enough to 
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establish that G > G(L + ) in [—1, 1]. Note that this inequality, as well 
as the notion of local minimizer, is independent of adding a constant 
to G. Hence, we may assume that 

G(s) = < G(L + ) for some s G [-1, 1], 

and we need to obtain a contradiction. Since G(L + ) > 0, we have that 

G(t) > e > for t in a neighborhood in [—1, 11 of L + 

1 + a w _ 

for some e > 0. 

Consider the points (b, 0, 0) = (xi — b, x^ — 0, . . . , x n — 0, y — 0) on 
dWf 1 . Since for R > 0, 

E B + (bm {u)> I -—G(u(x,0))dx 
fl -/r° (6,o) 1 + a 

and m(x,0) — L + , we deduce 

X\— )- + oo 

lim £ B + (6j0)0 )(u) > c(n) eiT for all R > 0. (5.2) 

6— >+oo 

The constant c(n) depends only on n. 

The lower bound f 15 . 2 jl will be a contradiction with an upper bound 
for the energy of u, that we obtain using the local minimality of u. 

For R > 1, let £r be a smooth function in M n+1 such that < £r < 1, 



1 in B7, ,„ 



on R n + +1 \B^ 
and |V£_r| < C{n){rjR)~ l , where 77 G (0, 1) is to be chosen later. Let 

y) := f + 6,ar 2 , ... , x n , y). 

Since 

(1 - ^R,b)U + £ft,&S = U + £ R , b (s - U) 

takes values in [—1, 1] and agrees with u on P^(&, 0, 0), we have that 

E B+ (6,0,0) ( U ) < ^B+(6,0,0)( M + ^,b( S - 

Next, we bound by above this last energy. Since G{s) = 0, the 
potential energy is only nonzero in B^ \ Bt, R , which has measure 
bounded above by C(n)i]R n . On the other hand, since we proved in 
Lemma 4.8(i) of [S] that 

\\VxU\\ L °°(B+( x ,a)) -> as x i ->■ ±0 °> 



2(» 
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we deduce that 

1 2 



lim / y°|V{u + &, 6 (s-u)}| 2 <2 / y a \V£ 

^ +00 JB+ (6,0,0) JB+ 
Ctrl') f R cm+l+a 

£ «-* = c(») w = o( B) . 



^2 

Putting together the bounds for Dirichlet and potential energies, we 
conclude that 

lim E B + (bm {u) < lim E B + ib00) {u + £ R)b {s - u)) 

< C{r]R n + V - 2 R n - 2s }, 

for some constant C > depending only on n, a, and G. 

Recalling the lower bound (15.21) . we now choose rj small enough so 
that Cr] = (l/2)c(n)e. In this way, (15 .2p and the last upper bound lead 
to (l/2)c(n)eR n < Cr]~ 2 R n ~ 2s . This is a contradiction when R is large 
enough. □ 

Proof of Theorem 12.111 We proceed exactly as in the proof of Theo- 
rem 1.4 in [H], page 1708. 

To prove part a), for R > 1 we consider problem (15. ip in a half-ball. 
Lemma I4TT1 gives the existence of a minimizer w with — 1 < w < 1. Note 
that in the lemma one needs condition (14. 5p . But in the presence of a 
layer, we showed in Lemma 4.8(i) of [8] that one has /(— 1) = /(l) = 0. 

On the other hand, Lemma 15.11 states that the layer u is the unique 
solution of (15. ip . Thus, u = w in B^. This shows that u is a local 
minimizer. 

To prove part b), G'{— 1) = G"(l) = was shown in Lemma 4.8(i) 
of [8]. We have established the other relation, G > G(—l) = G(l) in 
[— 1, 1], in Proposition 15.21 above. □ 

6. MONOTONICITY AND ID SYMMETRY OF STABLE SOLUTIONS IN 

R 2 . Proof of Theorem EUJ 

To prove Theorem 12.121 we need two lemmas. The following one, 
applied with d(x) = —(1 + a)^ 1 f'(u(x,0)), establishes an alternative 
criterium for a solution u of (II .3p to be stable. 

Lemma 6.1. Let d be a bounded and Holder continuous function on 
dRl +1 . Then, 



[ 2/ a |V£| 2 + / d(x)e>0 (6.i; 

Jr" +1 JdR n+1 
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for every function £ G C 1 (IR™ +1 ) with compact support in IR™ +1 , if and 



only if there exists a Holder continuous function ip in such that 

ip > in Wl +1 , (p G H^.(Rl +1 ,y a ), and 

L aV = in 

^ + d(x)<p = ondWl +1 . {b ' 2) 
dv a v + 

Proof. First, assume the existence of a positive solution ip> of (I6.2p . 
as in the statement of the lemma. Let £ G C 1 (1R™ +1 ) have compact 

support in IR™ +1 . We multiply L a ip = by £ 2 /V, integrate by parts and 
use Cauchy-Schwarz inequality to obtain ( 16. II) . 

For the other implication, we follow [9] . Assume that (16. ip holds for 
every £ G C l {W r l +1 ) with compact support in W]_ +1 . For every i? > 0, 
let Xr be the infimum of the quadratic form 



Qr(0= [ 2/lV£| 2 + f d(x)e (6.3) 



among functions in the class Sr, defined by 

S R = G H\B+, y a ) : £ = on T+ and £ 2 = l} 
C H 0>a (B+) = G H\B + R ,y a ) : £ = on T+} . 

We recall that the space Hq a {B^) was already defined in (I4.4p . 

By our assumption, A^ > for every R. By definition it is clear that 
Xr is a nonincreasing function of R. Next, we show that Xr is indeed 
a decreasing function of R. As a consequence, we deduce that \r > 
for every R, and this will be important in the sequel. 

To show that Xr is decreasing in R, note first that since d is assumed 
to be a bounded function, the functional Qr is bounded below in the 
class Sr. For the same reason, any minimizing sequence (£*,) has (V£fc) 
uniformly bounded in L 2 (B R ~,y a ). Hence, by the compact inclusion 
Ho }Cl (Bft) CC L 2 (Fjj.) (already mentioned in the proof of Lemma [4.ip . 
we conclude that the infimum of Qr in Sr is achieved by a function 

4>R G Sr. 

Moreover, we may take <pR > 0, since \<p\ is a minimizer whenever <fr 
is a minimizer. Note that 4>r > is a solution, not identically zero, of 

L a <pR = in B% 

+ d(x)(f)R = X R (f)R on T° R 

(p R = o on r£. 
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It follows from the strong maximum principle that <pR > in B R . 

We can now easily prove that Xr is decreasing in R. Indeed, arguing 
by contradiction, assume that Ri < R 2 and A^ = Xr 2 . Multiply 
La^Ri — by 4>r 2 , integrate by parts, use the equalities satisfied by 
4>r 1 and 4>r 2 , and also the assumption = Xr 2 . We obtain 



dv a 



>Ri 



0. 



and this is a contradiction since, on , we have <$r 2 > and the 
derivative d(j)R 1 /dh ,a < 0. 

Next, using that Xr > we obtain 



Jb + r Jr° R Jr° R Jr R 



for all £ G H 0ta (Bft), where 5r is taken such that < 5r < XR/\\d\\L°°- 
From the last inequality, we deduce that 

/ y a \V£\ 2 + [ d(x)e>e R [ y a \V£\ 2 (6.4) 
Jb+ Jr R Jb+ 

for all f G H 0A (B+), for e R > given by e R = 1 - 1/(1 + 5 R ). 

It is now easy to prove that, for every constant cr > 0, there exists 
a solution ifR of 

>+ 



L a <-PR = ° in B R 

+ d(x)(p R = on (6.5) 



dv a 

}Pr = c r on T^. 

Indeed, rewriting this problem for the function i[)r = <pr — cr, we need 
to solve 

X^b = inB + 



dip 

-q-^ + d(x)ipR_ + CRd(x) = on T° R 



ip R = on r^. 

This problem can be solved by minimizing the functional 



/ + y|V£| 2 + ( \\d{x)e + c R d{x)i\ 
Jb+ Jr° R I J 



in the space H a (B R ). Note that the functional is bounded below and 
coercive, thanks to inequality (16 .4p . Finally, the compact inclusion 
H 0ta (B R ) CC L 2 (T° R ) gives the existence of a minimizer. 
Next, we claim that 

(p R > in B^. 
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Indeed, the negative part ip^ of ipR vanishes on T^. Using this, ( 16. 51) . 
and the definition (16. 3p of Qr, it is easy to verify that Qr(p>~r) = 0. 
By definition of the first eigenvalue Xr and the fact that Xr > 0, this 
implies that ip^ = 0, i.e., p>R > 0. Now, the Hopf's maximum principle 
(Corollary 4.12 of [8]) gives <pr > up to the boundary. 

Finally, we choose the constant cr > in (16. 5p to have ^(0, 0) = 1. 
Then, by the Harnack inequality in Lemma 4.9 of [5] applied to p>s 
with S > 4R, we deduce 

sup <ps < Cr for all S > 4R. 

Now that (<ps) is uniformly bounded in B R , we use (12. 2p in Lemma [2721 
to get a uniform C^(B^ 2 ) bound for the sequence. Note that the 
constant C R in ( 12 .2p depends on the L°° (and not on the C a ) of dips, 
which we already controlled. However, to apply Lemma 12.21 we need 
to know that dps is C a . This is a consequence of the linear problem 
solved by ps an d the fact that dp>s € L°°. This leads to p>s G C a as 
shown in the beginning of the proof of Lemma 4.5 of [8]. 

Now, the uniform C^^B^^ bound gives that a subsequence of (ips) 

converges locally in to a C^ c (M™ +1 ) solution <p > of flO}. □ 

The previous lemma provides a direct proof of the fact that every 
layer solution u of ( II .ip is stable, which was already known by the 
local minimality property established in section 5. Indeed, we simply 
note that <p = u Xl is strictly positive and solves the linearized problem 
()672]), with d(x) = -(1 + a)- l f (u(x,0)). Hence, the stability of u 
follows from Lemma 16.11 

We use now the previous lemma to establish a result that leads eas- 
ily to the monotonicity and the ID symmetry of stable solutions in 
dimensions n = 1 and n — 2, respectively. 

Lemma 6.2. Assume thatn < 2 and thatu is a bounded stable solution 
of ( II. 3p . Then, there exists a Holder continuous function ip > in 
1R" +1 such that, for every i — 1, . . . , n, 

u Xi = Ciip in K™ +1 

for some constant Ci. 

Proof. Since u is assumed to be a stable solution, then ( 16.11) holds with 
d(x) := -(l + a)" 1 /' (u(x,0)). Note that d G C 13 by LemmaO Hence, 

by Lemma EH], there exists a Holder continuous function ip > in R™ +1 
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such that 

'L a ip = ® inM r ; +1 
(1 + a)~ 1 f\u{x, Q))(p = on dM$ 



dv a 

For % — 1, . . . , n fixed, consider the function 

u x 
a = — . 

<P 

The goal is to prove that a is constant in 
Note first that 

cp 2 Va = ipVu Xi - u Xt Vip. 

Thus, we have that 

div (y > 2 Vcr) = in M r | +1 . 
Moreover, we have 4^- = on cffi^ 1 since 

' av a t 

<p 2 (Ty = ifUy Xi ~ U X XPy = 

due to the fact that u Xi and ip both satisfy the same linearized boundary 
condition. 

We can use the Liouville property that we established in [8] (The- 
orem 4.10 of [8]), and deduce that o is constant, provided that the 
growth condition 

/ y a (^) 2 < CR 2 for all R > 1 (6.6) 

holds for some constant C independent of R. But note that (pa = u Xi , 
and therefore 



y\wf < / y a \^u 



Thus, we need to estimate this last quantity. 

To do this, we perform a simple energy estimate. Multiply the equa- 
tion div (y a 'Vu) = by £ 2 w and integrate in B^ Rl where < £ < 1 is 
a C°° cutoff function with compact support in B 2 r such that £ = 1 in 
B R and |Vf | < 2/R. We obtain 



y a {e\Vu\ 2 + 2^uV^Vu}= / (l + a)-V( M )eV 

Thus, by Cauchy-Schwarz inequality and since u and £ are bounded, 
/ y a e\Vu\ 2 < 1 - f y«i 2 \Vu\ 2 + C ! y a \Vi\ 2 + C\Tl R \ 
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for a constant C independent of R. Absorbing the first term on the 
left hand side, using that £ = 1 in B R and |V£| < 2/R, and computing 
fy R y a dy, we deduce 

/ y a | Vw| 2 < C{R- 2 R n R 1+a + R n } = C{R n ~ 2s + R n } < CR 2 
Jb+ 

since n < 2. This establishes (16.61) and finishes the proof. □ 
We can now give the 

Proof of Theorem 12.121 Let n — 2. The extension u of v is a bounded 
stable solution of (11.31) with / replaced by (1 + a)d~ 1 f. 

Lemma 16.21 establishes that u Xi = c^p for some constants q, for 
i — 1,2. If Ci = C2 = 0, then w is constant. Otherwise we have that 
— citt X2 = and we conclude that u depends only on y and on 
the variable parallel to (0, ci, c 2 ). That is, 

u(xi, rr 2 , y) = «o (cixj + c 2 x 2 )/{c\ + c 2 ) 172 , y) = u (z, y), 

where z denotes the variable parallel to (0,ci,C2). We have that u is 
a solution of the same nonlinear problem now for n — 1 thanks to the 
extension characterization; recall that the constant d s in (11.41) does not 
depend on the dimension. 

In particular 8 x uq = (c 2 + c 2 ,) 1 / 2 ^, and hence d x u > everywhere. 
This finishes the proof of the theorem. □ 

7. Layer solutions in R 

This section is devoted to the case n = 1. The Modica estimate that 
we proved in [5] (see Theorems 11.11 and 11.21 above) gave that 

G > G(-l) = G(l) in (-1,1) 

is a necessary condition for the existence of a layer solution in R. Note 
the strict inequality in G > G(±l). 

The rest of the section is dedicated to prove the existence of a layer 
solution under the above condition on G, in addition to G'(—l) = 
G'(l) = 0, as stated in Theorem 12. 41 The existence part of Theorem 12.41 
is all contained in the following lemma. 

Lemma 7.1. Assume that n = 1, and that 

G'(-l) = G'{1) = and G>G{-l) = G{\) in (-1,1). 

Then, for every R > 0, there exists a function ur G C^^B^) for some 
j3 G (0, 1) independent of R, such that 

— \<u R <\ in Bt, 



32 XAVIER CABRE AND YANNICK SIRE 

Mo,o) = o, 

d x u R > in B%, 
and Ur is a minimizer of the energy in B R , in the sense that 

E B +{u R ) < E B +(u R + ip) 

for every ijj G C 1 (B R ) with compact support in B R U T° R and such that 
-1 < u R + ip < 1 in B+. 

Moreover, as a consequence of the previous statements, we will de- 
duce that a subsequence of (u R ) converges in Ci oc (M.\) to a layer solu- 
tion u of (II. 3p . 

Proof. For R > 1, let 

Q+ = (-R,R) x (0,i? 1/8 ). 

Consider the function 

v (x, y) = v n {x) = — for (x, y) e QJ. 

arctan R 

Note that -1 < v R < 1 in Q%. 

Let u R be an absolute minimizer of E n + in the set of functions v G 
H 1 (Q R ], y a ) such that \v\ < 1 in and v = v R on d + Q\ in the weak 
sense. Since we are assuming G'(— 1) = = 0, the existence of 

such minimizer was proved in Lemma 14.11 We have that u R is a weak 
solution of 

'L a u R = inQ+ 

[l + a)^ = f(u R ) ond°Q + R 

u R = v R on d + Q R \, 

and, by the strong maximum principle and the Hopf 's lemma (Corol- 
lary 4.12 of 0), 

\u R \ < 1 in Q+. 

The function u R is Holder continuous by Lemma 12.21 

We follow the method developed in [9] and proceed in three steps. 
First we show: 

Claim 1: E Q +(u R ) < CR l/A (7.1) 

for some constant C independent of R. Here we take G — G(—l) = 
G — G(l) as boundary energy potential. We will use this energy bound 
to prove in a second step that, for R large enough, 

Claim 2: \{u\, 0) > 1/2}| > R 3/A and \{u%, 0) < -l/2}\ > R 3/4 . 

(7.2) 
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Finally, in a third step independent of the two previous ones, we prove 
that 

Claim 3: u R = d x u R > in Q+. (7.3) 

With the above three claims, we can easily finish the proof of the 
lemma, as follows. Since ^(-,0) is nondecreasing (here, this is a key 
point) and continuous in (—R,R), we deduce from (17.21) that for R 
large enough, 

u r (xr, 0) = for some xr such that \xr\ < R — R 3 ^ 4 . 

Since \x R \<R- R 3/4 < R - R 1/8 , we have that 

B+^ a (x R ,0) C (-R,R) x [0,i? 1/8 ] C Q+. 

We slide u R and define 

u R i/s(x,y) = u R (x + x R ,y) for (i,y)efl+ /1 (0,0). 

Then, relabeling the index by setting S = R 1 ^ 8 , we have that us £ 
C7^(5+(0,0)), -1 < u s < 1 in 5+(0,0), u s (0, 0) = 0, and d x u s > 
in Bg(0,0). Moreover, u,s is a minimizer in Bg(0,0) in the sense of 
Lemma 17.11 This follows from extending a given H l function if) with 
compact support in (Bg UT s )(xr, 0), and with |w+^>| < 1 in Bg(x R , 0), 
by zero in \ B^(x R ,0). Hence if) is a H l {Q^) function. Then one 
uses the minimality of u R in and the fact that the energies of u R 
and u R + ip coincide in \ Bg(x R , 0) to deduce the desired relation 
between the energies in Bg(x R , 0). 

Now we prove the last statement of the lemma: a subsequence of 
(u R ) converges to a layer solution. Note that we use the sequence (ur) 
just constructed, and not the sequence (u R ) in the beginning of the 
proof. 

Let S > 0. Since \ur\ < 1, Lemma 12721 gives C l3 (Bg) estimates for 
ur, uniform for R > 2S. Hence, for a subsequence (that we still denote 
by ur), we have that ur converges locally uniformly as R — > oo to some 
function u G C^ oc (M.\). By the additional bound (12. 3 j) on y a u y given 
by Lemma 12.2} one can pass to the limit in the weak formulation and 
u weakly solves (11 .3p 

We also have that \u\ < 1, 

m(0,0) = and u x > in R 2 + . 

Since u(0,0) = 0, we have \u\ ^ 1 and hence \u\ < 1 in , by 
the strong maximum principle and Hopf's lemma. Note that ±1 are 
solutions of the problem since, by hypothesis, G'(±l) = /(±1) = 0. 

Let us now show that u is a local minimizer relative to perturbations 
in [—1,1]. Indeed, let S > and if) be a C 1 function with compact 
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support in £><t U T° s and such that \u + < 1 in B$ . Extend ip to 
be identically zero outside B$, so that ip G H^ffi^). Note that, since 

— 1 < u < 1 and —1 < u + ip < 1, we have —1 < u + (1 —s)ip < 1 in B$ 
for every < e < 1. Hence, by the local convergence of (ur) towards 
u, for R large enough we have B$ C -B^ and — 1 < «r + (1 — e)^ < 1 
in S^, and hence also in B B . Then, since ur is a minimizer in 5^, we 
have E b +{ur) < E b +{ur + (1 — e)^) for i? large. Since -0 has support 

R R 

in Bt, U Tg, this is equivalent to 

E B + (u R ) < E B + (u R + (1 - e)if>) for R large. 

Letting i? — > oo, we deduce E B +{u) < E B +(u + (1 — e)^). We conclude 
now by letting e — > 0. 

Finally, since m x > 0, the limits L 1 * 1 = lim^^-i-oo 0) exist. To 
establish that u is a layer solution, it remains only to prove that iJ 1 = 
±1. For this, note that we can apply Proposition 15.21 to u, a local 
minimizer relative to perturbations in [—1, 1], and deduce that 

G > G{L~) = G(L + ) in [-1,1]. 

Since in addition G > G(— 1) = G(l) in (—1, 1) by hypothesis, we infer 
that l-Z^I = 1. But w(0,0) = and thus u cannot be identically 1 or 
— 1. We conclude that L~ = — 1 and L + = 1, and therefore u is a layer 
solution. 

We now go back to the functions u R defined in the beginning of the 
proof, and proceed to establish the three claims made above. 

Step 1. Here we prove (17.11) for some constant C independent of R. 
We take G — G(— 1) = G — G(l) as boundary energy potential. 

Since Eq+(u r ) < Eq+(v r ), we simply need to bound the energy of 

v R . We have 

|V^| = \d x v R \ = 1 —— — *— < C- ' 



arctan R 1 + x 2 1 + x 2 



and hence 



since < 1 + a < 2. 

Next, since G e C 2 ^, G'(-l) = G'{1) = and G(-l) = G(l), we 
have that 

G(s) - G(l) < C(l + cos(tts)) for all s G [-1, 1], 
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for some constant C > 0. Therefore, using that tt/ arctan i? > 2, we 
have 

G(v R (x, 0)) - Gil) < C { 1 + cos ( 7r arctanx ) 

I V arctan RJ 

< C(l + cos (2 arctan x)) 

2C 

= C2 cos 2 (arctan x) 



1 + x 2 
We conclude that 

\ {G(v R (x,0))-G(l)}dx<C -r— 2 <C. 

J-R J-R 1 + x 

This, together with the above bound for the Dirichlet energy, proves 

m- 

Step 2. Here we prove (17.21) for R large enough. 

Since u R = v R on {y = i? 1 / 8 } and J^ R v R (x) dx = 0, we have 



/H pH pH p 

u R (x,0)dx= u R (x,0)dx- u R (x, R 1/s ) dx = - I 
-R J-R J-r Jo: 



u R 



The energy bound ( 17. ip and the hypothesis that G — G(l) > give 
that the Dirichlet energy alone also satisfies the bound in (17.11) . We use 
this together with the previous equality and Cauchy-Schwarz inequality 
(writing \u R \ = y~ a / 2 y a / 2 \u R \), to deduce 



R 

u R (x, 0)dx 

R 



< I K\<\ / y' a - / y a ^ uRl2 

Jq + r I Jq+ Jq+ 

< C^RR^-^R 1 ^ 2 < CR 3 '\ 



1/2 

(7.4) 



since < 1 — a < 2. 

Next, by flU} we know that f* R {G(u R {x, 0)) -G(l)} dx < CR 1 ^ < 
CR 3 / 4 . On the other hand, G(s) - G(l) > e > if s e [-1/2, 1/2], for 
some £ > independent of R. Moreover, G — G(l) > in (—1, 1). We 
deduce 

, ri B/ /~\ \ i _ - l n l I / 

e\ v 

-R 



r R 

\{\u R (-,0)\ < 1/2}| < / {G(u R (x,0))-G(l)}dx<CR 3/ \ 

J-R 



and therefore |{|m*(-,0)| < l/2}\ < CR 3/4 . This combined with (17T4D 
leads to 

/ u R {x,0)dx <CR 3/ \ (7.5) 

y(-Ji^)n{|u«(-,0)|>l/2} 
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We claim that 

\{u R (; 0) > 1/2} | > R 3/4 for R large enough. 
Suppose not. Then, using (17.51) and \{u R (-,0) > 1/2}| < i? 3 / 4 , we 



obtain 
1 



\{u R (-,0) < -1/2}| < / u K (x,0)cte 



2' 



R( 



, < Ci? 3 / 4 . 

'(-fl,fl)n{te- B (-,0)<-l/2} 

Hence, all the three sets {|w R (-,0)| < 1/2}, {w R (-,0) > 1/2}, and 
{u R (-,0) < —1/2} would have length smaller than Ci? 3//4 . This is a 
contradiction for R large, since these sets fill (—R,R). 

Step 3. Here we establish the monotonicity result (17. 3p . This is done 
exactly as in Step 3 in the proof in [9], to which we refer. One simply 
uses the sliding method with the aid of the Hopf boundary lemma of 

0- □ 

Proof of Theorem 12.41 The necessary conditions on G follow from our 
previous paper [8]; see Theorem 11.21 above. 

That the conditions are sufficient for the existence of a layer v = v(x) 
follows from Lemma ITTj which gives a layer solution u = u(x, y) of the 
corresponding nonlinear extension problem (II. 3p and then by taking 
v := u(-,0). Note that we consider the extension problem with / 
replaced by (l + a)^" 1 / due to the relation (jl.4p between the fractional 
Laplacian and the Neumann derivative. 

Finally, the proof of the uniqueness result follows exactly that of 
Lemma 5.2 in [9] for the half-Laplacian. It uses the sliding method 
combined with the maximum principle Lemma 4.13 and Remark 4.14 
in our previous paper [8]. □ 

Proof of Theorem U .91 The proof is identical to that of Proposition 6.1 
in [9], page 1727. □ 
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